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ABSTRACT 


It is widely published that the terahertz (THz) spectral range has potential for 
imaging in the fields of military and security applications. The Sensors Research 
Laboratory previously achieved real-time imaging of concealed objects using a ImW 
quantum cascade laser (QCL) and an uncooled vanadium oxide/silicon nitride based 
microbolometer.. This thesis introduces an amorphous silicon based microbolometer 
with improved NETD in the 8-12 pm infrared spectral range. The QCL is usually 
operated in pulsed mode with rate in the hundreds of kHz which is much higher than the 
cut-off frequency of microbolometers of about tens of Hz. This indicates that neither 
camera should be able to detect the individual pulses of the THz beam. A detailed 
analysis showed that microbolometers can only detect the average power. Earlier 
experiments were then reproduced using the amorphous silicon based camera to assess 
the image quality but found it to be inferior to the silicon nitride based camera. These 
observations indicate that the absorption of THz in amorphous silicon is much weaker 
than silicon nitride. Other materials used to conceal military assets were analyzed and 
imaged to prove in principal the possibility of active THz imaging detection at a distance 
in narrow atmospheric windows. 
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I. INTRODUCTION 


In the days following the September 11, 2001 attacks on the United States the 
world security and threat environments changed drastically. Paralleling the political 
response was a push for technological solutions to address vulnerabilities that had been 
exploited by the attackers. Solutions included technologies leveraging our existing usage 
of the electromagnetic spectrum for security screening. Few of these threats were 
effectively countered, and new technological developments continued to be sought. One 
technology being considered was terahertz (THz) radiation, 0.3 to 10 THz, and in 
particular, THz imaging as a means of detecting hidden objects concealed by clothing. In 
2003 practical employment of such a system was far from being realized, and early 
papers exploring the topic were limited to proof-of-principal experiments promoting the 
potential of using non ionizing radiation with submillimeter resolution to perform 
security screening [1]. The principals were sound, and research into the use of THz 
radiation for security purposes quickly grew to encompass possibilities for detecting 
weapons of mass destruction and mine field detection through a combination of far IR 
spectroscopy and THz imaging [2]-[4]. By 2005, research in this area of THz detection 
of weapons and explosives was still largely theoretical [5], but papers analyzing THz 
transmission through various materials like common clothing and optical components 
continued to support its probable realization [7]-[8]. Experimental research was slow to 
evolve largely to a lack of THz sources and effective detectors for imaging and not from 
limited potential [9]. One research area developing between 2004 and 2006 involved the 
use of a quantum cascade laser (QCL) as a THz source in a THz imaging arrangement 
with a micrololometer camera as the detector [10]-[15]. The Sensors Research 
Laboratory (SRL) at the Naval Postgraduate School began THz imaging experimentation 
in this timeframe and advanced the pioneering work done at Agilent Laboratories in 
Santa Clara. Since that time, the SRL has contributed four research papers to the field, 
and the research reported in this thesis builds on that foundation [15]-[18]. 
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II. BACKGROUND 


A. IMAGING USING TERAHERTZ SOURCES 

THz imaging using QCL's holds significant promise for security and screening 
applications, largely due to its high average power, non-ionizing radiation and sub-mm 
spatial resolution. Early imaging work [10]-[15] provided meager results, but THz 
technology has continued to attract attention with advances in photonics and 
nanotechnology [19]-[20]. QCL performance was once a limiting factor, but THz QCL's 
can now operate at room temperature [21] and at power levels in the tens of mW's [19]- 
[22], The Sensors Research Laboratory, headed by Professor Gamani Karunasiri, has 
made contributions to the field of THz imaging using uncooled microbolometer cameras 
as detectors and QCL's as sources. Some significant images and videos published in 
scientific literature have come from LCol Barry Behnken and Professor Karunasiri. 
These can be seen in two papers titled "Real-time imaging using a 2.8 THz quantum 
cascade laser and uncooled infrared microbolometer camera" [17] and "Optimization of a 
3.6-THz Quantum Cascade Laser for Real-Time Imaging with a Microbolometer Focal 
Plane Array" [18]. An example of these images and the imaging arrangement can be 
seen in Figure 1 [18] and Figure 2 [17]. 



Figure 1. THz imaging arrangement. THz QCL is housed in a cryostat with emerging 
divergent beam collimated with parabolic mirrors passing through an object, and 
then focused onto a microbolometer focal plane array (from [18]). 
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Figure 2. 2.8 THz imaging of utility knife wrapped in tape shown in (a), (b) Single 

frame image produced by IR-160 camera, (c) 50 frame average, (d) 50 frame 

composite image (from [17]). 

The images were produced by an optically modified microbolometer camera, the 
IR-160 by Infrared Solutions. The significance of this work is that a very low powered 
(ImW) QCL has enabled imaging of a knife blade concealed in tape, the closest any 
research has come to realizing the highly touted potential introduced earlier. 

Where that research ended is essentially where this thesis research began. In an 
attempt to improve the imaging quality of our experimental set-up, other contributions to 
the field were studied. Most notable was a paper titled "Real Time Imaging Using a 4.3- 
THz Quantum Cascade Laser and a 320x240 Microbolometer Focal Plane Array" written 
by Alan Lee at the Massachusetts Institute of Technology [12]. That research group 
improved their results by synchronously modulating the QCL and bolometer focal plane 
array, thereby allowing differential imaging and reducing 1/f noise. This thesis research 
initially set out to accomplish differential imaging employing our pre-existing set-up. To 
do this a newer, more responsive digital vice analog camera was incorporated in the 
apparatus. The original and newer cameras are described in detail below. 
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B. INFRARED CAMERAS 


The camera used in the previous experiments [15]-[18] was an I n fr ared Solutions 
IR-160 un-cooled microbolometer camera [23], The camera is pictured in Figure 3. 



Figure 3. Infrared Solutions, Inc. IR-160 thermal imaging camera (from [23]). 

The focal plane array (FPA) consists of 160x120 pixels with a pitch of 50pm 
each. The FPA is comprised of a vanadium oxide (VOx) fdm applied to a silicon nitride 
(Si 3 N 4 ) membrane. This active layer is attached via support legs to a Si 3 N 4 substrate that 
acts as a heat sink. The structure can be seen in Figure 4. 



Figure 4. Schematic of an individual pixel from Infrared Solutions, Inc. IR-160 thermal 

imaging camera (from[34]). 
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The thermal response time is 10 ms, and the temperature eoeffieient of resistanee 
(TCR) is -0.023 k '. The eamera has a 30 Hz frame rate and is designed to operate in the 
8 to 14 pm speetral range with a noise equivalent temperature difference (NETD) of <100 
mK. The optics used in the IR regime is a 20 mm focal length, f/0.8 germanium lens 
with antireflective coating. The modified optics consisted of a f/1 Tsurupica lens with no 
antireflective coating. 

The second camera, pictured in Figure 5, is a Miricle 11 OK uncooled 
microbolometer made by Thermoteknix Inc [24]. 



Figure 5. Thermoteknix Systems Ltd. Miricle 1 lOK model thermal imaging camera 

(from [24]). 

Its FPA consists of 384x288 pixels, each with a pitch of 35 pm. The active and 
substrate layers are amorphous silicon with a thermal response time of 7 ms and TCR of 
-0.02 °C'\ The camera operates at a frame rate of 50 to 60 Hz and is designed for the 
7-14 pm, spectral region with a NETD of <50 mK. The optics used in the IR regime is a 
25 mm, f/I germanium lens with antireflective coating. The modified optics consisted of 
a f/1 Tsurupica lens with no antireflective coating. The camera also offers an external 
sync mode that allows synchronization of the camera with an external signal. This 
feature was necessary for experiments utilizing a differential imaging scheme. 

C. EARLY EXPERIMENTATION AEFECTS RESEARCH DIRECTION 

Early experiments conducted by simply swapping out the original camera failed 
to produce any images at all, let alone improved images. This directly prevented us from 
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achieving our original goals and heavily influenced the direction the thesis research took. 
The arrangement had to be re-analyzed, and the source was considered first. Previous 
experiments assumed that the best operating frequency was 300 kHz and that impedances 
were correctly matched. Careful analysis showed that there was an impedance mismatch 
affecting the output. The procedure for determining the best operating parameters for 
thermal imaging were developed using a 3.6 THz QCL and 2.8 THz QCL. This was the 
first major experiment undertaken. The results are reported in the following chapters. 
This still failed to yield improved results, and, in fact, the results were marginal 
compared to any of the previous work. Hoping to find answers in the theory, a careful 
review revealed that the theory presented in earlier work predicted that we should not be 
able to see an image at all given the optical power of the QCL. A theory was developed 
that adequately explained why we should be able to see images as demonstrated in 
previous research, but it did not explain the poor image quality then being attained. The 
cameras were then compared directly with each other, with surprising results. In spite of 
the much lower NETD for the new camera, <50 mK vice <100 mK, it produced very 
faint images when compared to the original camera, even with an identical beam and, 
identical optical components. This led to the idea that, although the NETD predicted 
better performance in the 7 to 14 pm range, the focal planes were different materials and 
perhaps they were responsible for the unpredictable results. The absorption 
characteristics of the focal plane arrays were then analyzed at THz frequencies and found 
to predict the bolometer output. The development of a plausible theory explaining why 
we are able to image at THz frequencies and the discovery of the significantly different 
focal plane responses at THz frequencies from IR wavelengths became the major focus of 
this thesis research. 

D. RESEARCH OBJECTIVES 

The objectives of this thesis research are as follows: 

1. Present a theory that adequately explains why we are able to see an image 
illuminated by a THz QCL using an optically modified un-cooled microbolometer. 
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Failures in the existing theory will also be presented through a eareful analysis of thermal 
imaging deviee figures-of-merit, namely responsivity, noise equivalent power and noise 
equivalent temperature differenee. 

2. Develop a method for optimizing QCL operation for use in THz imaging. 
Report results for impedance matching and frequency selection of the applied signal to 
yield the highest possible optical power output. Using the optimized operating 
parameters, investigate and report the laser IV characteristics for a 2.8 and 3.6THz QCL. 

3. Investigate THz propagation through optical components, material shown to be 
of interest in previous research and material of military significance presented in this 
thesis. Report results of proof-of-principal experiment as to whether THz imaging holds 
any promise for battlefield sensing. 

4. Demonstrate the performance of the IR-160 and Miricle 11 OK infrared 
microbolometer cameras in a THz imaging role. Contrast the performance of the two 
cameras in the IR regime, and show how the FPA material is responsible for the altered 
performance in the THz regime. Show how the FPA material also accounts for altered 
image quality at different THz frequencies. Make recommendations for an improved 
THz sensor based on this. 
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III. THEORY - SOURCE, DETECTOR AND PROROGATION 

MEDIUM 


A. THZ SOURCE - QUANTUM CASCADE LASERS 

The idea behind quantum cascade lasers (QCL) was first proposed by Kazamov 
and Souris in 1971, published in a paper titled “Possibility of amplification of 
electromagnetic waves in a semiconductor with a superlattice,” [25]. Although QCLs are 
based on semiconductors, the theoretical laser contrasts sharply to that of semiconductor 
lasers. Semiconductor lasers rely on electron hole recombination or interband transitions 
to emit photons. In bipolar semiconductor devices the photon wavelength is determined 
by this irradiative recombination across a band gap in the active layer of a forward-biased 
pn junction [26]. In contrast to this, the QCL is a unipolar device that uses a quantum 
well (QW) structure to emit photons when an electron undergoes an intraband transition 
within the well. When a voltage is applied to a series of spaced QWs, the electron can be 
made to "cascade" from one QW to a lower potential QW by consecutive tunneling, 
emitting a photon with each cascade. Implementation of practical QCLs has allowed 
lasers to operate at frequencies that are impossible with bipolar junction lasers, due to the 
extremely small band gaps required, and complications arising from thermal 
repopulation. The first QCL was produced at Bell Labs in 1994 by J. Faist et al [26]. 
Since then, QCLs have become the dominant mid-infrared semiconductor laser 
technology. In 2001 a laser with energy less then optical phonon energy was 
demonstrated at 67 pm, or 4.4 THz [27]. Selection of QW width and spacing, and 
advances in molecular beam epitaxy (MBE) fabrication techniques allowed the QCL to 
have a lasing frequency tailored by design, which has allowed operation to extend into 
the THz regime [22]. Steady advances from 2001 to today include: continuous wave 
(CW) operation at THz frequencies [28], low threshold currents [29], and high optical 
power with minimal cooling at THz frequencies [30]-[31]. This has opened up a world of 
possibilities for applications in numerous fields and has directly contributed to the 
advancement of THz imaging. 
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1. Quantum Cascade Laser Theory 


The key feature that distinguishes the QCL from other semiconductor lasers is the 
intersubband transition. A simple QCL model is a three level laser, with states 1 3), 12 ), 
and 11 ). Figure 6 shows how electrons are injected into the level 1 3> via the miniband in 
the left supperlattice barrier. The electrons that are unable to tunnel into the continuum 
via the minigap on the right fall to the level I 2) by emitting a photon. 


h[K'ife> 



The population inversion required for gain is created by the lifetime T 32 being 
significantly longer then the lifetime 121 . To achieve this, the lifetime T 32 for a level 
|3)->-|2> transition is prolonged due to the significant momentum transfer that must 
occur, and the lifetime 121 for a level 12 > 1 ) transition is reduced by making it resonant 

with the optical phonon energy of the host semiconductor. The resulting power achieved 
is the same as that for a conventional laser but with the addition of N active regions in the 
multiple quantum well (MQW) structure as seen in equation 1.1 [32]. 

P..=n„N—(0.1) 

q 
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2.8 and 3.6 THz QCLs 


The lasers used in our imaging experiments are a bound-to-continuum design 
optimized for 2.8 and 3.6 THz operation. The epitaxial MQW layers consist of 
Alo. 15 Gao. 85 As and GaAs with 120 periods for the 3.6 THz laser and 120 for the 2.8 THz 
laser. The band diagram and layer widths (in nanometers) of the quantum wells are 
shown in Figure 7, with the GaAs layers highlighted in bold characters. 



D-di.'ped 


lli/3.2/11/2.4/ I L/1.S/I2/U / I VR/l ‘’1(V/0.9/]f..4.-0.t/V/A.5 

(a) (b) 

Figure 7. (a) Band structure and (b) waveguide of bound-to-continuum 3.6 THz QCL 

used in this experiment (from [ 10 ]). 

It is very evident from Figure 7 that, if the applied bias is not correct, it is very 
likely that electrons will not tunnel into their intended energy state and that the device 
will not function as intended. This will seriously deteriorate the performance 
characteristic we are most interested in, the optical power output. The device must be 
cooled to approximately 10 K to reduce thermal excitation of electrons, and it heats up 
significantly with the necessary applied currents of close to 1 A. Characterization of THz 
QCL performance for optical imaging purposes must consider all these factors. They are 
discussed in detail below. 

3. 2.8 and 3.6 THz QCL Performance 

There are a number of factors that contribute to the performance of a QCL. Band 
misalignment and power reduction can occur from an improperly applied bias, impedance 
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mismatch between power supply, and laser as well as thermal effects. Thermal effeets 
are mitigated through cryogenie cooling and pulsed mode operation, while the device's 
dynamie impedance can only be determined experimentally through 1-V measurement. 
The lasers used in this researeh were analyzed in previous works authored by NPS 
students. The first analysis coneentrated on the impedanee issues involved in optimizing 
power output [15], while the others analyzed effects related to both the pulse duty eyele 
and the operating temperature [16]-[18]. 


An assignment of unique impedance eharacteristie for the device is difficult due 
to the dynamic nature of its operation. It does not have a conventional "operating point". 
The 3.6 THz QCL impedanee was previously reported as 3.47 Q by Agilent Laboratories 
[10]. The SRL construeted a step-down transformer using this information to match the 
known 50 Q eharacteristie impedance of coaxial cables [15]. The transformer turns ratio 
appropriate to the deviee was ealeulated by eonsidering the transfer eharaeteristics. The 
power (Pin) to the primary windings of a transformer and the power (Pout) out of the 
seeondary windings must be matehed for optimal power transfer between source and 
load. Since the power can be expressed as P=IV or P=I^Z where I is eurrent, V is the 
voltage and Z is the impedance it follows that ZinIin^=Zout lout^- It should be noted that 
power losses are redueed when eurrent is reduced. The induced EMF (^tum) is the same 
in the primary and seeondary windings, allowing it to be expressed as ^tum=Vp/Np=Vs/Ns 
where Vg and Vp are the voltages in the seeondary and primary windings, and Ng and Np 
are the number of windings in the seeondary and primary windings. This can all be 
expressed in the following equation. 


7 =7 

^out 




vNgy 


( 0 . 2 ) 


The transformer eonstrueted used the turns ratio (Np/Ng)^ = 0.26, based on the 
laser's reported 3.47 Q. impedance. Follow-on work reported by the SRL [16]-[18] used 
this transformer and determined that the impedance was dynamic and depended on 
eryostat temperature and duty eyele, but, despite these variations, this transformer was 
ultimately used for both the 2.8 and 3.6 THz QCLs. 
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The duty cycle can have an effect on both the device temperature and the device 
impedance. This can be seen in Figure 8, where the I-V characteristics for this device 
were measured at a 300 kHz pulse repetition frequency and over 5%-25% duty cycles. 



(a) (b) 

Figure 8. (a) I-V analysis for the 3.6 THz QCL operated under 5, 10, 15, 20, and 25% 

duty cycles, (b) I-V analysis for the 3.6 THz QCL operated at 10% duty cycle 
for cryostat temperatures ranging from 0 K to 100 K (from [18]). 

The impedance was measured to be 2.7 Q [18], which differs substantially from 
that reported in Reference [10]. Finally, the cryostat temperature can affect the dynamic 
impedance. This can be seen in Figure 8. where the I-V characteristics for this device 
were measured at a 300 kHz pulse repetition frequency with a 10% duty cycle and over a 
temperature range of 0 K-lOO K. The dynamic impedance can clearly be seen to be 
changing again in this analysis. Optimized QCL performance for this laser clearly has to 
focus on narrowing the dynamic impedance to that corresponding to peak THz power 
output, the region where optimal imaging would occur. 

The performance optimization in reference [18] focused on analyzing the duty 
cycle for peak output power. Figure 9 shows how the duty cycle affects the peak output 
power as measured by a FTIR for duty cycles of 5, 10, 15, 20, and 25% . 
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Figure 9. (a) Peak power analysis for the 3.6 THz QCL, operated under 5, 10, 15, 20, 

and 25% duty cycles and (b) peak power analysis for the 3.6 THz QCL, operated 
under 10, 15, 20, and 25% duty cycles over time (from [18]). 


The peak power was delivered with a 20% duty cycle upon initial start up of the 
laser cooled to 10 K. The cryostat temperature, and therefore device temperature, was 
affected by the duty cycle. The temperature change was then shown to have an effect on 
the signal intensity. Over a 4 minute lasing interval the signal intensity would degrade. 
The best operating duty cycle maximized the power applied to the detector while 
minimizing Joule heating of the device. Figure 9 clearly shows this optimization 
occurring at 20%. A 15% duty cycle did not put as much power on the detector, but at 
25% heated up the device's temperature increased significantly through Joule heating, 
which degraded the output. 

In all of the above analyses, a 300 kHz pulse rate was employed, and the 
transformer described earlier was used. As it applies to this project, the results of the 
20% duty cycle output power optimization was consistently used, but all other 
assumptions about the optimal frequency, proper dynamic impedance range and 
transformer construction were not. These issues are addressed in the experimental 
procedure, section IV, and were a new method for determining the best operating 
parameters for both the 2.8 and 3.6 THz QCLs were developed. 
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B. 


DETECTOR - MICROBOLOMETER 


Uncooled microbolometer focal plane arrays (FPAs) for thermal imaging were 
first developed in 1983 by Honeywell using silicon micromachining and the well known 
dependence of material resistance on temperature (bolometric effect). The project was 
funded by both the US Army Night Vision Laboratories and the Defense Advanced 
Research Projects Agency. In 1989 they made a device consisting of 64 x 128 pixels 
consisting of silicon nitride (SiN) membranes and vanadium oxide (VOx) resistive film 
with a temperature coefficient of resistance (TCR) of approximately 2% per °C. The 
detector was attached to a monolithic read out integrated circuit [33]. Since then, 
uncooled solid state infrared imaging arrays have continued to advance with material and 
micromachining breakthroughs. Although microbolometers FPAs of this type were 
developed primarily for infrared imaging, the bolometric effect is responsive across the 
entire electromagnetic spectrum. These devices can therefore be used to detect radiation 
at all wavelengths and, in particular, have applications in terahertz imaging. 

The figures-of-merit typically used to quantify the performance of infrared 
microbolometers can be adopted to quantify the device performance in the terahertz 
spectrum as the devices are adapted for use in this regime. 

1. Microbolometer Theory 

As mentioned above, a microbolometer is a thermal detection device that uses 
heat flow to affect a measurable quantity in a sensing element. The radiant power Pq 
incident on the detector causes the sensing element to undergo a temperature increase 
proportional to the heat absorbed by the pixel at that particular position in the array. In 
the case of a microbolometer, this temperature change AT causes a change in resistance 
AR. If the absorbed radiation is small, this change is linear with the incident power [34]. 
The degree with which this effect is felt is quantified by the TCR, a, and the device 
resistance R. This is expressed as 

AR = aRAT. (0.3) 

Figure 10 shows schematic of a typical microbolometer pixel structure with a thin film 

vanadium oxide (VOx) resistive sensing element. 
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Figure 10. Schematic of a typical microbolometer pixel structure (from [34]). 

Proper thermal isolation of these structures prevents leakage conduction to 
neighboring pixels and convection losses to the surrounding atmosphere. In the ideal 
case, the temperature change on the sensing area is followed by heat flow down the 
support legs to the substrate. If the sensing area has a heat capacity C (J/K) the relation 
between the heat Q the element absorbs with a temperature change AT is Q = CAT. The 
change with respect to time is therefore dQ/dt = C(dAT/dt). As the temperature of the 
element rises, with respect to the substrate, heat conduction begins. The absorbed 
thermal energy is then lost through the support structure. This loss is characterized by a 
thermal conductance G (W/K), and it occurs at a rate GAT. When radiative power Po is 
incident on the sensing area, assuming that the device remains approximately in 
thermodynamic equilibrium, the heat balance equation can be expressed as 

1 A 

P„=GAT + C-. (0.4) 

dt 

In practice, it is desirable to modulate the incident radiation beam, in order to employ 
standard noise-minimization techniques such as synchronous detection (lock-in methods), 
dynamic averaging (boxcar integration), 1/f noise suppression, and so forth. Hence, Po is 
deliberately forced to be an explicit function of time, Po(t). The solution of equation 1.4 
may then be expressed as a time-shifted linear superposition of the detector's impulse 
response (Green's function), which is exponential with a characteristic time constant x 
[34]-[35]. This can be seen in Figure 11. 
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Time 

Figure 11. Modulated radiation ineident on foeal plane array. 


When modulated radiation power Pq e*®* is ineident on the sensing area with 
absorption efficieney r\, neglecting average power, the heat balance equation can now be 
expressed as 

qPoej“‘ = GAT + C ^. (0.5) 

By assuming a Green's function of the form e^^‘ , where the thermal time constant 
X = C/G, we can solve the above equation as follows: 


r Aatc^'Mi = ^ f e^''"^j“^‘dt, 

Jdt C J 


As time goes to infinity, or in the steady state, the result is 

qPoe^“‘ 


AT = Re 


(G + jcoC)e 


j'i> 


( 0 . 6 ) 


(0.7) 


where tan((p) = cox. This result can be recast as [35] 

AT =—sin((at-(|)). (0.8) 

gVi + co^" 

From this we can gain an idea of how a pixel should respond when operated within the 
context of pulse-modulated terahertz imaging. If we assume a pixel to be thermally 
isolated, with typical thermal parameters C = 1x10'^ J/K and G = 1x10’’ W/K, the thermal 
time constant is then x = 10 ms [36]. The peak power of the QCL used in previous 
experiments [15]-[18] is about 2.5 mW when operated at a pulse repetition frequency of 

300 kHz. The effective power was degraded by the use of a Tsurupica window and a 
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Tsurupica lens in the beam each with approximately 64% transparency at 2.8 THz. For 
the time being, we will only consider one attenuating element, keeping in mind that this 
preliminary discussion represents the absolute upper performance limit for our system. 
Under these conditions the pixel temperature change in a 19200 pixel FPA is calculated 
to be about 

AT = 44 pK. (0.9) 

It should be noted at this point that this is a liberal estimate, since the 20% operational 
duty cycle of the QCL was not factored in, and the power incident on and absorbed by the 
FPA was assumed to be 100% of the incident flux which is not the case. To gain an 
appreciation for how the cameras should nominally respond to such signals, they must be 
compared to the figure-of-merit noise equivalent power (NEP). The derivation of NEP is 
made in the following sections, but the result for the IR-160 camera is 230 pW. The 
associated temperature change per pixel necessary to produce an image pixel with unit 
signal-to-noise ratio is then NEP/G, which is AT = 2.3 mK. That implies that it is 
impossible to detect individual laser pulses at 300 kHz using a microbolometer pixel with 
a AT of only 44 pK. However, due to the relatively long time response x compared to the 
laser's pulse period, we clearly detect measurable signals from the laser due to a 
cumulative effect from multiple pulses. Therefore the analysis needs to include not only 
the high frequency modulation component of the laser power, but also the quasi-static 
contribution coming from the effective averaging produced by the detector's slow 
response. 

2. Analysis of Microbolometer Operation under Pulsed Power 

When equation 1.8 was derived it was assumed that the average power is small, 
which is obviously not correct in reality. It is simply a consequence of linear small-signal 
analysis, commonly used in all formal derivations requiring only the frequency response. 
A more complete representation of Figure 11 includes the quasi-static power, expressed 
as Po, and the time varying power, Poe”®* . The heat balance equation 1.4 can now be 
expressed as 


18 



( 0 . 10 ) 


P„ + P„ej“‘ =GAT + C-. 

“ “ dt 

This differential equation can be solved separately for each of the power terms (quasi¬ 
static and time varying) on the left hand side. Since we have already solved the 
modulated power term we will solve the unmodulated average power case, add the 
modulated solution and apply the boundary conditions. It follows that 


AT-Pq/G 


r 1 

J —dt and ln(AT - P^ / G) = -t / x + const. Solving for AT we are left 


with the following. 


AT = ^ + Be' 
G 


( 0 . 11 ) 


Again, not all power incident on the device will be thermally absorbed depending on the 
absorptive properties of the material in a particular frequency regime and an efficiency 
factor r\ must be applied. Adding this to the previously solved modulated radiation yields 
the amount of the temperature change in a pixel as. 


AT = + pe'^ + — , ' ° sin(03t - (|)). 

G GVl + coV 


( 0 . 12 ) 


Applying the boundary condition AT = 0 at t = 0 we find P = 


expression for the temperature change in a pixel is. 


riPp sin((^) qPp 

gJI^" g • 


AT = -^ 1-e ^ ^ sin((Dt -(|)) + e ^ sin (|) . (0.13) 

G V / Gvl + co^x^ I y 


As time goes to infinity we are left with the final form of equation 1.14. 

AT = ^ + —Ji^siii(o)t-(|)). (0.14) 

G Gv1-i-co^t^ 

The second term in equation 1.14 is negligibly small at typical pulsing frequencies of the 

QCL, and the microbolometer response is essentially coming from the quasi-static term 

due to average power of the laser. The second term response (solid line) can be seen in 

Figure 12 for Pq = 1.25 mW, r\ = 0.64, G = 10’’ W/K and x = 10 ms in a 19200 pixel 

array. Comparing this response to that of the first, average power term (dashed line) and 
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it is evident that beyond a pulse repetition frequency of 1 kHz the average power term 
dominates as the time varying term approaches 0. 



Figure 12. Temperature change in a pixel for Po = 1.25 mW, r\ = 0.64, G = 10'^ W/K and 
T = 10 ms. Response of second term in Equation 1.14 (solid line) compared to 
average power term (dashed line). 

In order to gain an immediate understanding of the solution presented above we 
have used a technique to model the thermal and electrical properties of a microbolometer 
pixel in PSPICE [37]. The final form of the heat balance equation for the above solution 
is given as, 

1 A'T' 

qPf, + riPoej“‘ = GAT + C-. (0.15) 

dt 

Using a simple RC circuit given in Figure 13 (a), the thermal parameter G (W/K) 
equivalent is the inverse resistance, R (1/G), and the thermal parameter C (J/K) 
equivalent in the circuit is the capacitance, Cc (F). The power equivalent is the AC 
current I(t) (A), and the temperature is therefore taken as the voltage response V(t) (V). 
The heat balance equation is now presented in terms of electrical equivalents as 

l(t) + I(t)eJ“‘ =-iAV(t) + C,^^. (0.16) 

R dt 

The thermal time constant x = C/G is now represented as an RC time constant x = RCc. 
Using the thermal parameters C = 1x10'^ J/K, G = 1x10 ’ W/K and the thermal time 
constant x = 10 ms [36] we calculate the electrical equivalents as R = 10 Mf2, Cc = 1 nF 

and X = 10 ms. The peak power of the QCL used in previous experiments [15]-[18] is 
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about 2.5 mW when operated at 300 kHz. As calculated earlier, the associated 
temperature change per pixel in a 19200 pixel array was AT = 0.42 K. This value was 
taken as the applied voltage in the circuit, with the equivalent voltage offset. Next the 
circuit was behavior was modeled using analog modeling in PSPICE. Figure 13. (b) 

shows the temperature rise and decay, respectively, when the laser is on and off at a 
modulation frequency of 300 Hz for several time constants of a microbolometer pixel. 
At this frequency the time varying term is evident in the pixel response. Figure 13 (c) 
depicts the circuit response at 300 kHz. At such a high frequency, the rapid modulation 
response dies out, leaving only the quasi-static response. It shows, that for standard 
video rate, at 30 frames per second, THz imaging can be done using pulsed QCLs and 
micobolometers, provided that the average power of the laser is appreciable. 



(a) 



Figure 13. (a) RC circuit used to model equation 1.16. (b) 300 Hz PRF and (b) 300 kHz 

PRF. 


21 










3. 


Responsivity St 


The responsivity St, is defined as the ratio of a sensors output signal to the input 
stimulus, whieh in this case is incident power. In the case of microbolometers, the change 
in resistance due to the temperature change is given by AR = aRAT. Substituting for AT 
using equation 1.13, the change in resistance can be obtained as 

sin((Bt -(])) + e sin (|) 


AR = 


aRqPf, 


G 


1-e 

V y 


+ 


aRriPg 


gVi + co"t" 


(0.17) 


The resistance change can be measured in a conventional Wheatstone bridge, as in 
Reference [36]-[37] using one microbolometer and three equal resistors. The current i 
through the side of the bridge with the device under an applied bias Vb is 




1 = 


and Ro(l + aAT) is the resistance through the device. Solving for 


Rq (1 + ctAT) + Rq 
output signal of the bridge, Vout, gives the following expression. 


RoVb 


V _ 

2 Ro(l + aAT) + Ro 


(0.18) 


Assuming aAT«l, then AV^^j, = (aVgAT)/4. Using this result, the responsivity is 


found by dividing Vout by the incident power, Pq, yielding 

7 .r r 


91 = 


aVgri 


4G 


1-e 


V 


+ 


aVgq 


4gVi + co"t" 


sin(cot - (|)) + e ^ sin (() 


(0.19) 


Using Vb = 5 V, a = 0.02 K'\ q = 0.64, and G = 1 x 10'^ V/W and all other parameters 
unchanged gives 

91 = 1.6x10'V/W. (0.20) 

A value of 91 = 2.0x10' V/W was previously reported in Reference [17] for the 8 to 12 
pm spectral region, using only the frequency response term for a similar camera. These 
responsivities would allow us to see an image using the average power of the laser. If an 
entire focal plane were uniformly illuminated with a 1.25 mW average power beam, and 
the FPA consisted of 19200 pixels, the fraction of the power incident on a single pixel 
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would be about 65 nW, ignoring losses between pixels. The resulting voltage output 
would be about 10 mV, whieh can be easily measured. 


4. Noise Equivalent Temperature Difference (NETD) 


The noise equivalent power (NEP) of a device is the hypothetical input radiation 
power Po that would be required to create a detected signal-to-noise ratio of 1. The noise 
in a resistive bolometer is dominated by Johnson noise if the readout bias is less then 100 
pA [37]. If this is not the case, other noise sources can be significant and may include 
thermal, Johnson, 1/f and amplifier noise. Previous work [16]-[18] calculated the NEP 
for the IR-160 microbolometer and used the result to develop a NETD, equation 1.21, for 
the THz spectral region. That earlier treatment assumed only Johnson noise and made 
non-applicable assumptions about the bolometer resistance. Because of this, and the fact 
that 1/f noise in amorphous silicon can vary by three orders of magnitude with growth 
technique [37], this method could not be applied to our Miricle camera. Instead, NEP 
was deduced for both cameras in the IR region and used to determine the NETD in the 
THz region. 


NETD is defined as the temperature change of the target that produces a FPA 
signal equal to the total noise [34]. The NETD can therefore be expressed as 

4F2 


NETD = 


• NEP 


( 0 . 21 ) 




T.A(dP/dT)^ 

where F is the f-number of the lens, A is the pixel area, and xo is the product of 
transmittance of the lens, absorptivity of the FPA, and the thermal emissivity of the target 
object. The differential spectral radiance dP/dT is the change in power per unit area 
radiated by a blackbody at temperature T, with respect to T, between some predetermined 
wavelength limits A,i and X 2 .. This differential is calculated from the spectral radiance, as 
given by Plank's Radiation Law, equation 1.22, and depicted in Figure 14. 

27rh v^ 


P(T) = - 


C2 ehv/kT_j 


( 0 . 22 ) 
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Figure 14. Spectral Radiance of a blackbody at 300K (from[16]). 


Typical values for dP/dT in the 3-5 pm and 8-14 pm range for a 295 K blackbody 
are 2.10x10'^ Wcm'^K’' and 2.62x10 '^ Wcm'^K'*, respectively [34], To determine this 
quantity in the terahertz range P(T) can be found from 

27ih v^ 


P(T)=f 

j\ 


7 .^hv/kX 1 

1 c e -1 


dv, 


(0.23) 


if it is assumed that the source is an ideal blackbody. Because we are at the low 
frequency limit of the blackbody spectrum, the Plank formula (1.23) may well be 
approximated by the Rayleigh-Jeans formula. Since hv/kT«l and x= hv/kT is small, 
e’^-l+x, and we are left with the following: 


(•V' 

P(T)=[ 

Jv, 


27i:h 


v 


c 1 -I- hv / kT -1 


dv. 


(0.24) 


P(T)= f'^^v^dv. (0.25) 

Jv, c 

To a good approximation, the lower limit of integration can be set to 0. Setting the upper 
cut-off frequency to be Vc, the integral in equation 1.25 is then reduced to 
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3c 


2 c 


The quantity dP/dT is given simply as 


dP 2nk 


dT 3c 


-V 


2 c • 


(0.26) 


(0.27) 


This value is found to be 4x10'^ Wcm'^ K'*, using Vc = 5 THz. Substituting the result into 
equation 1.21 then gives the following expression for NETD at typical THz frequencies 
[16]-[18] : 


6c^F^ 

NETD =-^ • NEP. 

ToATlkv, 


(0.28) 


Previous analysis of the IR-160 camera [16]-[18] yielded an NETD of 3 K with a 
cut off-frequency of 5 THz. This contrasts significantly with the NETD of 0.1 K that is 
typical in the IR spectral region. The increase in NEDT indicated the need for additional 
illumination for imaging in the THz frequency range. 

As mentioned earlier, the NETD of 3 K in the THz region was based on NEP 
calculations made with approximations. Since the NETD is known for both cameras in 
the IR spectral region along with dP/dT of 2.62x10“^ Wcm'^ K''[34] in that region, we 
can make a direct calculation of the NEP which remains constant in all regions. The 
NETD of the IR-160 camera is <100 mK with a f/0.8 germanium lens and 50x50 pm^ 
pixel. The target emissivity is taken as 1, the germanium lens with anti reflective coating 
has a transmittance of 90%, and the SiN mean absorptivity is 80% [38]. Because of the 
fact that the camera was engineered for almost complete IR absorption in the 7 - 14 pm 
range, we will take its absorptivity as I, but this will not necessarily be the case for the 
THz regime. Using these values in equation 1.22 gives a NEP in the 7 to 14 pm region of 
230 pW for the IR-160 camera. The Miricle camera has a 35x35 pm pixel, a NETD of 
<50 mK and a, f/1 germanium lens with anti reflective coating. Its NEP is calculated to 
be 36 pW. 

Using these values an accurate NETD for the THz spectral region can be 
extrapolated. The Tsurupica lens used for THz imaging is a f/1 untreated lens with a 
mean transmission of 0.60 in 1 to 5 THz range. Because the microbolometers used were 

25 



not designed for operation in the THz region the absorptive properties of the FPA eannot 
be assumed to be unity. The mean absorptivity in the THz region for SiN (used in the IR- 
160) is 40% and for Si (used in the Miriele 11 OK) is 10%. These values, along with a 
dP/dT of 4x10'^ Wcm'^ K ' reported earlier, give NETDs of 38 K for the IR-160 eamera 
and 49 K for the Miriele camera. This predicts that the Miriele camera should be better at 
imaging in the IR spectral region with diminished performance in the THz regime 
compared to the lR-160 camera. In the next section, terahertz propagation through the 
atmosphere and different military relevant materials are presented. 

C. THZ TRANSMISSION IN THE ATMOSPHERE AND MATERIALS 

Losses in uniformly absorbing material can be characterized by Beer's Law 

P = P„e““, (0.29) 

where a is the absorption coefficient per unit length x of propagation through a material. 
At THz frequencies the effect can vary significantly between materials and frequency 
across different intervals. In a THz imaging system, carefiil operating frequency 
selection is critical for standoff imaging. A delicate balance must be struck between 
undesirable attenuation, from optical components and the intervening atmosphere, and 
desirable attenuation created by the target material. Each of these components will be 
addressed separately. 

1. Atmospheric Transmission 

Probably the most discouraging effect that limits the feasibility of stand-off THz 
imaging is the atmosphere, where water vapor is the dominant absorber that attenuates 
THz power drastically [39]. Several computer codes like LOWTRAN, MODTRAN and 
HITRAN have been developed to quantify atmospheric effects [40]. For propagation 
modeling of narrow bandwidth radiation, high resolution is required and for this 
FASCODE has been developed. FASCODE accesses the HITRAN database directly and 
many of the models in LOWTRAN indirectly like, molecular continuum absorption, 
molecular scattering, and aerosol absorption and scattering but at a much higher 
resolution [40]. Figure 15. was generated from FASCODE with a resolution of 100 
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MHz at an altitude of 100 m on a clear day with 50 km visibility. The path was aim 
horizontal column at sea level and the frequency band was between 1.5 and 5 THz. 



1.5 2.0 2.5 3.0 3.5 4 0 4 5 5.0 

Frequency [THz] 


Figure 15. FASCODE model of atmospheric loss at sea level on a clear day. 

For imaging done in this thesis two lasers of 2.8 and 3.6 THz were used. 
Originally being considered was a demonstration of stand-off imaging, but Figure 16. 
shows why this was not feasible due to strong atmospheric absorption at these 
frequencies. Imaging distances had to be carefully limited in all experimentation. 
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(a) (b) 

Figure 16. FASCODE model of atmospheric loss for a) 2.8 and b) 3.6 THz windows. 

The FASCODE model in the THz region may not be entirely correct. The 
HITRAN model by itself gives fairly accurate predictions compared to FTIR analysis of 
the atmosphere as shown in Reference [14]. In Figure 17 (a) we can see that the losses 

27 






predicted by HITRAN at 296 K and 40% relative humidity are significantly less then 
those predicted in Figure 15 using FASCODE. The regions expanded in (b) and (c) show 
the results of FTIR analysis (dashed line) under the same conditions compared to the 
HITRAN model (solid line). Figure 18 (a) and (b) are provided as a comparison to the 
FASCODE model. 
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Figure 17. a) HITRAN model and FTIR results of atmospheric loss for b) 4.3 and c) 4.9 
THz windows. HITRAN (solid line) and FTIR results (dashed line). Lasing 
spectra of QCL's are depicted with dashed vertical lines (from [14]). 
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Figure 18. FASCODE model of atmospheric loss for a) 4.3 and b) 4.9 THz windows. 
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THz Laser Beam Propagation 


To gain an appreciation of the feasibility of stand-off detection using a 
QCL/microbolometer imaging combination, the propagation of the laser beam through 
the atmosphere has to be considered in conjunction with the responsivity of the camera. 
Later we will see that a 1.25 mW average power beam is sufficient to image through an 
18% transmissive target material and two 64% transmissive optical components. The 
images produced were at the lower limit of our experimental sensitivity, so it can be said 
that the lowest power incident on the FPA to allow imaging with the IR-160 camera is 92 
pW. Before we consider beam size and intensity, we can use this lower limit of 
detectable power to determine the absolute upper range limits possible for stand off- 
imaging. Due to the extreme atmospheric attenuation at the operating frequencies of our 
current QCLs, propagation distances will be estimated for the next laser to be used by 
the SRL, a 200 mW, 4.9 THz QCL. Using equation 1.29, the attenuation is calculated 
over distance for the FASCODE model, HITRAN model, and measured FTIR extinction 
coefficients shown in Figure 19. 



Figure 19. Maximum propagation distance for 200 mW 4.9 THz QCL using FTIR 
measured, FASCODE modeled and HITRAN modeled extinction coefficients. 

When the beam divergence and losses due to optics are considered, the distances 
shown in Figure 19 to permit detection of a 92 pW beam at the FPA will, of course, be 
reduced. Because of the structure of QCLs (see Figure 7 (b)) and their long operating 

wavelengths relative to the sizes of optical components, the beam emitted can be 

29 







divergent due to diffraction. It is necessary to collect and then to collimate the beam with 
parabolic mirrors to obtain an approximately collimated beam capable of travelling any 
significant distance. 

3. FPAs, Optical Components (Tsurupica) and Plastic Tape 

The next most critical component in our THz imaging system was the 
microbolometer FPA. Early in the experimentation process it was believed that the lower 
NETD of the Miricle camera in the infrared range would yield better results, but it was 
quickly apparent that poor absorption efficiency of THz radiation by the FPA pixels 
generated inferior images. Figure 20 (a) and (b) show the transmittance of THz through 
SiN (used in IR-160 camera) and amorphous Si (used in Miricle camera), respectively. It 
can also be noted that the transmittance differs markedly with frequency. The IR-160 
FPA is a VOx film on SiN, and it is evident in Figure 20 that at both 2.8 and 3.6 THz this 
material is drastically more absorptive than the amorphous Si used in the Miricle 11 OK 
camera. Interestingly, the figure also explains the improved imaging attained by the IR- 

160 camera at 3.6 THz from that of 2.8 THz in previous work done by the SRL [16]-[18]. 
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(a) SiN and (b) Si (after [41]). 

The next materials analyzed were the Tsurupica used for the cryostat windows 
and camera lenses and the plastic tape used to obscure the knife blade in previous 
experiments [16]-[18]. The Tsurupica absorption was important to quantify since there 

were two components in the imaging system, the window and lens made of this material. 
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Figure 20. Percent absorption of FPAs. 






Materials were analyzed using a Nicolet Magna 850 Fourier Transform Infrared (FTIR) 
spectrometer. The spectrometer generates FIR radiation from an Ever-Glo blackbody 
light source in the 700 to 50 cm'' wavenumber range. The internal beam splitter was a 
solid substrate, and the detector was a deutereated triglycine sulfate (DTGS) pyroelectric 
detector. Samples were introduced into the internal sample holder and scanned at a 
mirror velocity of 0.6329 cm/s for a total of 32 background scans and 32 scans of the 
sample. The resolution of the scans was 4 cm '. It was not experimentally viable to scan 
the focal plane material of the two cameras, SiN and amorphous Si elements. Data 
previously obtained from a similar FTIR spectrometer was accessed from the RIKEN 
THz database [41]. 


Figure 21 (a) and (b) show measured transmittances using our FTIR for Tsurupica 
and Polytape, respectively, in the 1.5 to 6 THz range. Only about 60-64% of the beam 
power (Po) passes through these components at 2.8 and 3.6 THz frequencies. For this 
reason, optical components like lenses had to be used sparingly, and highly 
THz-reflective gold plated paraboloidal mirrors had to be used to steer and to manipulate 
the THz beam. Also seen in Figure 21 (b) is the highly THz-transparent quality of the 
tape used to conceal the knife blade in previous measurements [16]-[18]. A significant 
amount of success was achieved imaging metallic objects through this material. 



Figure 21. Percent transmission of (a) Tsurupica 4.0 mm, and (b) polytape 0.25 mm. 
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To further advance this work new and more challenging materials were sought. 
Of particular interest was a class of materials designed to conceal military assets that 
were specifically designed to defeat detection systems operating in other areas of the EM 
spectrum. These materials are discussed below. 

4. Camouflage and Suppression of Military Assets 

Some commonly used methods for concealing military assets are screens, 
blankets, heat reduction systems, and mobile camouflage nets. There is a significant 
amount of theory behind the operation of all these measures [42] that are designed to 
suppress detection across most areas of the electromagnetic spectrum, from UV to radar. 
Due largely to a lack of sources and detectors, the THz spectrum has not yet been 
exploited, nor have adequate countermeasures been developed. Some properties of THz 
radiation, like its ability to penetrate materials, make it difficult to counter with the 
systems mentioned above. When considered in conjunction with the inherently high 
reflectivity of metals and ceramics at THz frequencies, this leads to tantalizing 
possibilities for real-time imaging of concealed military assets from modest stand-off 
distances. A number of these systems were analyzed at THz frequencies in the results 
section, 5-B. 
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IV. EXPERIMENTAL SET UP AND THZ IMAGING 


A. OPTIMIZATION OF QCL POWER 


The first step taken towards studying improved THz imaging capability was 
enhancing the optical power emitted by the 2.8 and 3.6 THz QCLs. As mentioned 
earlier, previous analysis of laser performance concentrated on optimizing optical power 
by analyzing the duty cycle. The research reported in this thesis focused on optimizing 
the optical power by analyzing the dynamic impedance characteristics of the laser and 
matching impedances for the optimized operating frequency. The bracketing procedure 
outlined below yielded the optimal operating conditions for terahertz imaging. 


The laser itself is housed in a Janis closed cycle refrigeration (CCR) system 
vacuum shroud. The atmosphere is removed from the CCR shroud using an Edwards 
Picodry pump. The cold finger within the shroud is cooled using a high pressure helium 
CCR system and is maintained at 10 K with a LakeShore 312 auto-tuning temperature 
controller. To dissipate the heat generated by the large currents pumped through the 
device, the QCL was mounted on a copper block heat sink attached to the cold finger of 
the refrigeration unit, as illustrated in Figure 22. 



Figure 22. 3.6 THz QCL mounted on copper thermal dissipation block attached to cold 

finger. Red arrow in inset depicts direction of laser output beam. Shroud has 
been removed from the CCR system (from [16]). 
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The window of the shroud was a 4 mm thick Tsurupica plate whose THz 
transmission properties can be seen in Figure 21. The beam was collected and collimated 
by a 50.8 mm focal length gold plated off axis-parabolic (OAP) mirror, and steered with 
a gold plated flat mirror into the external port of the FTIR spectrometer. The QCL was 
driven with an Agilent 8114 100 V, 2 A pulse generator, equipped with a variable output 
impedance selection capability. The frequency was varied from 1 to 400 kHz, with 
impedances of 5Q, lOQ, 15Q, and 20Q. Impedances were further bracketed using a 
frequency that produced the most optical power as measured by the FTIR. Once the 
impedance matching of the signal generator to the load of the QCL was found, we then 
investigated whether the transformer used in previous experiments [16]-[18] provided a 
better impedance match and improved optical power output. This defined the impedance 
that would be used subsequently for improved imaging. 

The pulse repetition frequency (PRF) was then refined using the aforementioned 
optimal impedance value. Optical power measurements were taken for 1, 100, 200, and 
300 kHz PRF by varying the current from 0.5 A to 1.6 A for the 3.6 THz laser and 0.5 A 
to 1.4 A for the 2.8 THz laser. This identified the optimal PRF. 

Next the duty cycle was investigated briefly. Previous imaging using this laser 
determined that a 20% duty cycle was optimal for short duration lasing. This conclusion 
needed to be confirmed for our new optimal operating parameters. Previous studies also 
concluded that operating time affected performance. The effect on optical power was 
analyzed over 1 min and 5 min times using the new impedance parameters and the 
various frequencies investigated before. These results were plotted to determine the best 
operating frequency for short term and sustained imaging applications. 

The laser I-V characteristics were then taken for both the 2.8 and 3.6 THz lasers, 
based on their optimal settings found from the above procedure. The voltage and current 
were measured with an Agilent Infmiium Oscilloscope, and the optical output was 
measured concurrently from the FTIR spectrometer. 


34 



B. OPTIMIZED IMAGING ARRANGEMENT 


Imaging was achieved using the optical set up previously established by the SRL 
[16]-[18]. The optical arrangement can be seen in Figure 1. This thesis research 
optimized the optical path to achieve improved imaging. The exact procedure is 
described below. Both QCLs were driven with an Agilent 8114 pulse generator at 300 
kHz PRF and a 20% duty cycle. For the 2.8 THz QCL, a 1.3 A current was applied and 
impedance matched at 5 Q. The 3.6 THz QCL had an applied current of 1.4 A and was 
impedance matched at 13 f2. The QCLs were cooled to 10 K using the Janis CCR system 
described above. The shroud window was a 4 mm Tsurupica flat with transmission of 
64% and 60% for the 2.8 and 3.6 THz QCLs, respectively. The devices were placed as 
close to the window as possible within the shroud, and the highly divergent THz beam 
was captured using a 50.8 mm focal length gold plated OAP mirror. This mirror was 
placed inverted and close to the device outside of the shroud so an un-collimated beam 
emerged. The beam was sent to a second OAP m irror of focal length 101.6 mm. The 
second mirror was placed approximately 160 mm above the first mounted on the same 
optical post. A Tsurupica 1/f lens with a focal point of 20 mm was then attached to the 
camera to focus the image on the FPA as illustrated in Figure 23. The use of an of-axis 
paraboloidal mirror to collect the output radiation from the QCL and then concentrate it 
on the sample was not an optically optimized way to proceed. The better choice would 
have been an ellipsoidal mirror, but one was not available. 



Figure 23. Schematics of the optics used in THz imaging. 
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V. RESULTS 


A. OPTIMIZATION OF QCL PERFORMANCE 

The 3.6 THz QCL was driven by the pulse generator with 1.4 A, 10% duty cycle 
at frequencies from 1 to 400 kHz and impedances of 5, 10, 15, and 20 Q. Figure 24 (a) 
shows that peak optical power production is achieved somewhere between 1 and 100 kHz 
PRF and around 15 Q impedance. A frequency of 100 kHz was chosen to further bracket 
the impedance match that produced the most optical power as measured by the FTIR 
spectrometer. Figure 24 (b) shows this occurring at a 13 impedance. 



Figure 24. (a) Peak power of the device given a pulse generator impedance 50 Q into 5, 

10, 15 and 20 Q at various frequencies, and (b) peak power of the device at 100 
kHz given a pulse generator impedance of 50 into 11 to 15 Q . 

Knowing the optimal impedance match for the device, we verified the 
performance by again looking at the optical output for 1, 100, 200, and 300 kHz PRF. 
This was done for both the optimal 13 Q impedance and for the outlying impedances, 
since it is possible the device impedance may be frequency dependent. Figure 25 (a) 
confirmed that the optimal impedance match and frequency combination was 13 Q and 
lOOkHz PRF. This also yielded the optimal current as 1.55 A. Since previous 
experiments used a frequency of 300 kHz and a transformer to impedance match the 
device and driver, the new results had to be compared to prior performance, given these 


37 





different parameters. Figure 25 (b) shows a drastic improvement at a current of 1.55 A, 
100 KHz PRF and a 13 Q impedance. It also shows that previous experiments were 
achieving promising results with only 70% of the available optical power. 
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Figure 25. (a) Peak power of the device given a 10% duty cycle, and a pulse generator 

impedance of 50 Q into 13f2 for various currents and pulse frequencies, and (b) 
peak power of the device given a 10% duty cycle, and a pulse generator 
impedance of 50 into 13 Q or with transformer for various currents and PRF 

of 100 and 300 kHz. 

The 20% duty cycle was assumed to be optimal from previous experiments [16], 
but the optimization above had to be done at 10% to avoid damaging the device. The 
power available at 20% the duty cycle was investigated briefly as shown in Figure 26 
(a), where over a two-fold increase in optical power was obtained with operation at 20% 
vice a 10% duty cycle. Previous studies also concluded that there is an effect on 
performance with operating time. Figure 26 (b) shows the effect on optical power over 
time (1 min and 5 min) using a 13 Q impedance and the various frequencies investigated 
before. These results indicate the best operating frequency for short term and sustained 
imaging is again 100 kHz. 
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Figure 26. (a) Peak power of the device at 100 kHz for 10 and 20% duty cycles, with a 

pulse generator impedance of 50 into 13 Q at various currents, and (b) peak 
power of the device given a 10% duty cycle, with a pulse generator impedance of 
50 Q into 13 at various frequencies after 1 and 5 minutes of lasing. 

The above analysis was undertaken for the 2.8 THz laser, with minor 
modifications. Care had to be taken, and the device performance was not pushed as hard 
as possible, since at the conclusion of the analysis above the 3.6 THz laser ceased to 
operate and was no longer available for imaging experiments. Figure 27 shows the 
results of the iterative process that again yielded 100 kHz as the optimal operating PRF 


with a device impedance of 5 



Figure 27. (a) 2.8 THz QCL peak power given a 10% duty cycle, and a pulse generator 

impedance of 50 into 3 to 7 Q, and (b) peak power given a 10% duty cycle, 
and a pulse generator impedance of 50 into 5 Q for various currents and PRFs. 
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Finally the I-V characteristics of the two devices were investigated using the 
optimized parameters generated above. Figure 28 (a) and (b) show the I-V curves for the 
3.6 and 2.8 THz lasers, the threshold current for lasing and the applied current and bias 
that resulted in optimal optical power output. 




(a) (b) 

Figure 28. a) I-V characteristics and peak power output of the 3.6 THz QCL given a 

10% duty cycle, with a pulse generator impedance of 50 Q into 13 at a PRP of 
100 kHz, and b) I-V characteristics and peak power output of the 2.8 THz QCL 
given a 10% duty cycle, with a pulse generator impedance of 50 Q into 5 at a 

PRF of 100 kHz. 

B. ANALYSIS OF CAMOUFLAGE MATERIALS 

Combat clothing is the most prevalent material used to conceal personnel. The 
analyzed sample of combat uniform is comprised of two types, one for temperate 
woodland regions and the other for arid regions. Both samples are 0.32 mm thick and are 
pictured in Figure 29 (a). The transmission spectra in Figure 29 (b) shows little promise 
for THz screening of military personnel clothed in this material, due to poor THz 
transmission. 
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(a) (b) 

Figure 29. a) Photographs of combat uniform, arid region (top), 0.32 mm, and temperate 
woodland (bottom), 0.32 mm, and b) percent transmission of arid region combat 

uniform. 


The next sample is a garnish for a temperate woodland camouflage screen. The 
garnish mimicking foliage consists of a green and black patches each 0.49 mm thick. 
The sample is pictured in Figure 30 (a). Since these materials typically do not conceal 
large metallic or ceramic objects, they likely could not be defeated by THz imaging 
system due to poor THz transmission. 
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Figure 30. a) Photograph of temperate woodland camouflage screen, black and green 
garnish, 0.49 mm, and b) percent transmission of temperate woodland 

camouflage screen. 

The next sample is from a blanket vice a garnish and is intended for winter 
operations. The sample is pictured in Figure 31 (a). As with the woodland screen, it 
holds little promise for THz imaging due to poor THz transmission (see Figure 31 (b)). 
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Figure 31. a) Photograph of winter operations camouflage screen, 1.12 mm, and b) 
percent transmission of winter operations camouflage screen. 

Thermal blankets and heat reduction systems are commonly used to conceal metal 
and ceramic assets. The thermal blanket samples in Figure 32 (a) are for temperate 
woodland, arid regions and winter operations. THz transmission spectrum in Figure 32 is 
similar in all three samples and attenuates too drastically for practical imaging. 
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Figure 32. (a) Photographs of thermal blankets, temperate woodland (top left), arid (top 

right) and winter operations (bottom center), and (b) typical percent transmission 

for all samples. 




The Heat Transfer Reduction System is comprised of an insulator, a 0.44 mm 
base fabric, and 0.16 mm scrim. The arid region system without the base layer is 
pictured in Figure 33 (a). Figure 33 (b) shows that the heat transfer reduction system 
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holds promise for imaging below 2 THz if it is used without the insulating layer. When 
used with the insulating layer THz energy is totally attenuated. 
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Figure 33. a) Photograph of heat transfer reduction system without insulating layer, 
base layer (right) and scrim (left), and b) percent transmission of heat transfer 
reduction system, base layer and scrim. 

Camouflage nets are the most widely used means of concealing vehicles and other 
large assets. The analyzed sample is a garnish for temperate woodland camouflage 
netting used for vehicles. The garnish mimics foliage and consists of a green and black 
leaf pattern and a brown and black leaf pattern, each about 0.25 mm thick. The sample is 
pictured in Figure 34 (a). This material holds significant possibilities for THz imaging 
based on the transmittance spectra shown in Figure 34 (b). It possesses relatively good 
THz transmission properties and effectively conceals metallic and ceramic objects that 
reflect THz radiation strongly. The possibility of achieving reflective THz images from a 
system based on a UAV platform is very high. This material will be used for imaging 
studies in the experimental portion of the thesis. 
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Figure 34. 



(a) 


(b) 


a) Photograph of woodland vehicle camouflage net, and b) percent 
transmission of garnish. 


The next sample is a base layer and scrim for an arid camouflage system again 
used for vehicles. The material consists of a 0.50 mm base layer and 0.27 mm scrim. 
The sample is pictured in Figure 35 (a). Like the woodland vehicle camouflage net this 
material holds significant possibilities for THz imaging if the scrim is used alone. With 
the base layer employed it would defeat a THz imaging system. 
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Figure 35. a) Photograph of the arid region mobile camouflage system, base layer (left) 
and scrim (right) and b) percent transmission of scrim only. 




The last sample is a winter operations mobile camouflage system for vehicles. 
The system is comprised of a 3.43 mm backing, a 0.47 mm base layer and a 0.51 mm to 
0.54 mm garnish layer. The system is pictured in Figure 36 (a). When elements are 
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employed either individually or together as a system they will defeat THz imaging 
systems, based on the transmission spectra shown in Figure 36 (b). 



Frequency |THz] 



(a) (b) 

Figure 36. a) Picture of winter operations mobile camouflage system, garnish, backing 
layer, base layer, and garnish, and b) percent transmission base layer only. 

C. PERFORMANCE OF TWO CAMERAS AND THZ IMAGING 

Before performing any THz imaging, the respective performances of the two 
cameras were compared in the IR spectral region. Figure 37 (a) shows an optical image 
of the imaging setup used, along with the same view in the IR using the Miricle 11 OK 
(Figure 37 (b)) and IR-160 (Figure 37 (c)) cameras with their original optics. It can be 
clearly seen that the Miricle produces sharper image compared to the IR-160, due to its 
smaller pixel size and larger FPA format. 



(a) (b) (c) 

Figure 37. a) Photograph of the imaging setup (from [16]), b) IR image of (a) using 
Miricle 1 lOK camera, and b) IR image of (a) using IR-160 camera. 
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Next, the two cameras were compared at THz frequencies using this optical setup 
and the 2.8 THz QCL. The beam was focused directly on the FPA using OAPs without 
the Tsurupica lens attached to the camera. A current of 1.3 A was applied to the QCL at 
100 kHz PRF, 20% duty cycle and impedance matched at 5 Q. Figure 38 shows what 
was initially unexpected, the lower NETD camera (in the IR range) produced the 
weakest THz image. The THz absorptive properties of amorphous Si were compared to 
that of SiN (see Figure 20), and it was immediately obvious that we would have little or 
no success with the new camera due to weaker absorption efficiency of THz radiation by 
amorphous Si. Also apparent are the diffraction patterns above the beam which were the 
first indications that the imaging arrangement used in previous experiments [16]-[18] 
may not have had the optical properties originally thought. 



(a) (b) (c) (d) 


Figure 38. a) Background using Infrared Solution IR-160 camera, b) 2.8 THz beam 
using 50 Q into 5 at a PRF of 100 kHz, 20% Duty Cycle, c) background 
using Thermoteknix Miricle 1 lOK camera and d) 2.7 THz beam using 50 Q into 

5Q at PRF 100 kHz. 

The above findings brought the fact that FPA's behave very differently at THz 
frequencies than in the IR spectral region to the forefront of our investigations. As a 
second look at this important result, we looked at earlier published work done by the 
SRL. 

The images created using 2.8 and 3.6 THz frequencies of the knife blade 
concealed in tape pictured in Figure 39 show a frequency dependence of image quality. 
The 2.8 THz image was slightly weaker than the 3.6 THz image taken with the IR-160 
SiN based camera. From Figure 20, and assuming no reflection, the FPA should exhibit 
85% absorption at 3.6 THz compared to 70% at 2.8 THz. Figure 20 (a) and (b) also show 
the absorption characteristics at 2.8 THz for Si at 40% and SiN at 70%. Clearly, the 
amorphous silicon based camera was not a good choice for THz imaging, as seen in 
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Figure 40. The 2.8 THz beam is seen with the Tsurupica lens in place in (a), and the 
unconcealed knife (b), and concealed knife (c), are just barley visible in the bottom left 
hand comer of the beam. As a consequence, no further imaging was attempted with the 
Miricle camera. 



(a) (b) (c) 

Figure 39. a) Knife concealed by plastic tape, b) 2.8 THz image using Inf rared Solution 
IR-160 infrared microbolometer camera and c) 3.6 THz image using Infrared 
Solution IR-160 infrared microbolometer camera . 



(a) (b) (c) 


Figure 40. a) 2.8 THz beam, b) Un-obscured knife, and c) Knife obscured with tape. 

Next, the identical imaging scenario was attempted with the IR-160 camera. 
Figure 41 (a) shows an identical beam pattern to that of the faint Miricle beam pattern 
above. The pattern is similar to that of previous imaging at 2.8 THz shown in Figure 40 
(b). Interestingly, the results of optimizing the QCL can bee seen in these images. 
Figure 41 (b) shows the un-obscured blade, and the hole in the blade can clearly be seen, 
where that was not the case previously. Figure 41 (c) shows the blade obscured by 
plastic tape where unlike the image in 39 (b) the only indication that the blade was 
obscured in the new image is a slight change in the diffraction pattern most noticeable in 
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the lowest ring. In fact, the presence of such a diffraction pattern implies that some 
fundamental assumptions about the optical setup may be questionable. This is discussed 
below, in the conclusions. 



As a final experiment, the possibility of using this technology in a battlefield 
sensing environment was investigated. After analyzing numerous materials used for 
concealing military assets, the camouflage net was considered as possibly possessing 
THz vulnerabilities. At 2.8 THz the material was found to be 18% transmissive. The 
optical arrangement used to produce the above images was again used, and a small 
camouflage leaf was placed on the blade. Figure 42 shows how the knife is clearly 
visible under the leaf 



(a) (b) (c) 

Figure 42. Knife covered with woodland camouflage net a) camouflage net coming into 
beam, b) covered knife coming into beam, c) covered knife fully visible in beam. 


The significance of this image cannot be overstated. The ability to image through 
a material only 18% transmissive, at an average power of 1.25 mW, demonstrates the 
incredible potential of this technology. This proof-of-principle experiment clearly 
indicates the possibility of active THz imaging detection at a distance using an 

arrangement designed to operate in narrow atmospheric windows. 
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VI. CONCLUSIONS 


The main aim of this thesis research has been to investigate how infrared 
optimized microbolometer absorbing membranes interact with THz frequencies and 
whether it will be feasible to perform active stand-off detection of concealed targets using 
THz radiation. In addition, a set of commonly used camouflage materials were also 
probed to determined their effectiveness and vulnerabilities in the THz spectral range. 

During the course of research, it was found that in THz imaging using relatively 
slow microbolometer sensors and pulsed QCL sources, the most important parameter is 
the average power output of the source. The QCL is usually operated in pulsed mode 
with rate in the hundreds of kHz which is much higher than the cut-off frequency of 
microbolometers, which is typically of about tens of Hz. This indicated that neither 
camera should be able to detect the individual pulses of the THz beam. A detailed 
expression was derived to understand the response of microbolometer pixels to pulsed 
QCL sources, which clearly shows that the pixels can only respond to average power. 

This directly influenced the need to reconsider QCL operation. Armed with this 
knowledge, a method was developed to optimize peak optical power by maximizing the 
average power that would be sensed by the microbolometer. The average power was 
optimal for both the 2.8 and 3.6 THz QCL's at 100 kHz PRF, and with impedances 
matched at 5 Q and 13 Q, respectively. 

A direct comparison of the optically modified IR-160 and Miricle 11 OK 
microbolometer cameras yielded surprising results. Earlier experiments were reproduced 
using the amorphous silicon based Miricle 11 OK camera to assess the image quality, but 
we found it to be inferior to the silicon nitride-based IR-160 camera. These observations 
indicated that the absorption of THz in amorphous silicon is much weaker than in silicon 
nitride. This was confirmed with an FTIR analysis of the two FPA materials. The 
differing camera response is predicted based on the absorption efficiency of the two 
materials. 
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In addition, for stand-off imaging applications, it is very important to select THz 
frequencies that propagate well through the atmosphere with minimal extinction. Correct 
frequency selection allows for operation in narrow atmospheric windows where 
attenuation is relatively low and a high average power fluency can be maintained 
between the camera and target. A frequency that minimizes attenuation through optical 
components and concealing materials, but interacts strongly with concealed objects must 
also be chosen correctly. Finally, based on this frequency selection, a focal plane array 
material can be engineered to achieve maximum absorption efficiency and, therefore, 
optimal imaging response. This was shown in the image quality differences at 2.8 and 
3.6 THz when imaging using the IR-160 silicon nitride-based microbolometer. As we 
have seen in this thesis, improved microbolometer performance in the IR spectral range 
does not necessarily translate into enhanced performance in the THz spectral range, and 
that there are technical limits to the success of modified IR microbolometers as THz 
imaging detectors. 

Finally of interest to the author was the notion that current camouflage 
technologies for concealing military assets might easily be exploited by a THz imaging 
system. Of the numerous materials analyzed only one showed promise, the camouflage 
net and it happens to be the mainstay for concealing mobile military assets by countries 
around the world. The clear images of the knife concealed by the THz transparent single 
camouflage leaf can easily be extrapolated to an artillery piece concealed by an entire 
camouflage net. The outstanding results from this proof-of-principle experiment indicate 
that THz imaging holds promise for battlefield sensing so long as improvements to 
sources and detectors continue to be made. 

There are several avenues for further research suggested by our work. In 
particular, it is evident from the images of the concealed knife blade that the experimental 
setup could be improved. In particular, the images exhibit pronounced background 
interference fringes and/or diffraction patterns, which seriously degrade the target field- 
to-background signal-to-noise properties of the images. If their interference effects can 
be understood and then suppressed during later experiments, significantly greater target 
discrimination should result. Several interferences may be drawn from the images. First, 
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notice in Figure 39, 41 and 42 narrow regions of the target knife blade do not appear to 
produce edge diffraction fringes. This suggests that the high-contrast fringes surrounding 
the blade image are artifacts of the optical setup, rather than the target object. Second, 
the fringes appear to have an eccentric geometry, rather than one that is circularly 
symmetrical about the beam axis. This, as well as the total absence of linear Haidinger 
fringes, indicates that optically flat components such as windows, are not the culprits. 
Third, the fringe eccentricity strongly resembles the phase-space footprint that would be 
produced by an off-axis conic section. Two such components exist in our setup, namely 
the paraboloidal mirrors. These mirrors are standard components, manufactured for use 
at much shorter wavelengths than those used here. At wavelengths on the order of 1 mm, 
the optical coatings on these m i r rors are extremely thin, much thinner than the radiation's 
skin depth. We speculate then, that the background interference fringes arise from the 
off-axis paraboloidal mirrors. By using alternative components machined from solid 
metallic blocks, such as bulk aluminum, it may be possible to suppress the fringes. This 
could be tested quite easily because, at millimeter wavelengths, high quality optical 
surfaces can be fabricated using standard machining methods. Expensive optically 
figured surfaces are not needed at all. Indeed, such machined components could replace 
the 90° off-axis configuration shown in Figure 23, which is far from ideal, because the 
paraboloids' focal surfaces are severely inclined to the optical axis. The testing of these 
hypotheses could form the basis for follow-on research to this thesis project. 
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